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We apply the Jaynes principle of maximum entropy for the partial reconstruction of correlated 
spin states. We determine the minimum set of observables which are necessary for the complete 
reconstruction of the most correlated states of systems composed of spins-1/2 (e.g., the Bell and the 
■ Greenberger-Horne-Zeilinger states). We investigate to what extent an incomplete measurement 

' can reveal nonclassical features of correlated spin states. 
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I. INTRODUCTION 



The seminal paper by Vogel and Risken JÏJ on the tomographic reconstruction of Wigner functions of light fields has 
greatly enhanced interest into the old problem of "measurement" of states of quantum-mechanical systems. Within 
the last few years tomographic reconstruction has been experimentally realized for example by Raymer and coworkers 



and Mlynek and coworkers Q. Tomographic reconstruction schemes of states of other bosonic systems such 
I ■ as vibrational modes of trapped àtoms Q and atòmic waves [[3] have been proposed. Recently Wigner functions 
of vibration states of a trapped atom have been experimentally determined by Wineland and coworkers |>|, while 
f*»» ■ Kurtsiefer and coworkers have measured Wigner functions of atòmic wave packets. Leonhardt || has extended the 
0^ ' ideas of Vogel and Risken to the case of Wigner functions in discrete phase spaces associated with physical systems 
, with finite-dimensional Hilbert spaces, such as spin systems. 

The problem of reconstruction of states of finite-dimensional systems is closely related to various aspeets of quantum 
information processing, such as reading of registers of quantum computers This problem also emerges when states 
of àtoms are reconstructed. In particular, Walser, Cirac, and Zoller jïo| have shown that under certain conditions 
d | quantum states of a single quantized cavity mode can be completely transferred on to the infernal Zeeman submanifold 
of an atom. Consequently, the reconstruction of the states of a cavity mode is reduced to the problem of reconstruction 
of angular momentum states in a finite dimensional Hilbert space. 

From the postulates of quantum mechanics it follows that the complete reconstruction of a state of a quantum- 
mechanical system can be performed providing a complete set of system observables (i.e., the quòrum [|ïl"||) is measured 
on the ensemble of identically prepared systems. This goal on onc hand may be technically difficult to realize and on 
the other hand it may not be necessary. In many situations even partial knowledge (i.e., incomplete reconstruction) 
of the state is sumeient for particular purposes. 

The complete reconstruction of spin states have been addressed in the literature |]8|, |Ï2| , p]JÏ3| . In the present paper 
we will analyze the problem of the partial reconstruction of these states. We will show how the spin state can be 
reconstructed when just a restricted set of mean vàlues of the system observables is known from the measurement. 
Utilizing the Jaynes principle of maximum entropy we will partially reconstruct the density operator from the available 
(i.e., measured) mean vàlues of system observables. 

We will also address the question as to which is the minimum observation level (i.e., a specific subset of system 
observables) on which the complete reconstruction can be performed. In particular, we will analyze the reconstruction 
of the most correlated states of the system composed of two and three spins-1/2 (i.e., the Bell and the Greenberger- 
Horne-Zeilinger states) . 

The present paper is organized as follows. In Section II we briefry review the principle of maximum entropy and 
the formalism of state reconstruction associated with particular observation levels. In Section III we present a simple 
illustration of a reconstruction of a state of a single spin- 1/2. Section IV will be devoted to the detailed analysis of 
the state reconstruction for a system of two correlated spins-1/2. In Section V we will address the problem of the 
(partial) reconstruction of the Greenberger-Horne-Zeilinger states. In Appendix A we present detailed reconstruction 
of density operators on two nontrivial observation levels. 
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II. RECONSTRUCTION OF DENSITY OPERATORS OF QUANTUM STATES 



When it is a priori known that experimental data contain the complete information about the state of the system, 
then it is just a question of technical convenience how to perform a transformation of this data into a more familiar 
object suen a density operator. A particular example of this procedure is quantum homodyne tomography JÏJ when 
from the measured probability distributions of rotated quadratures one can reconstrucij^] (with the help of the inverse 
Radon transformation) the Wigner function of the state. 

Now we can ask the question: "What is the density operator of the quantum mechanical system when an incomplete 
measurement over this system is performed? " . In this case the experimental data does not provide us with sufficient 
information to specify the density operator of the system uniquely, i.e. there can be many density operators which 
fulfi.il the constraints imposed by incomplete experimental data. In this situation one can only estimate what is the 
most probable density operator which describes the system. 

In principle we can distinguish two different forms of incompleteness: firstly, when a precise knowledge of a subset 
of the quòrum JÏT[ of system observables is known; secondly, when system observables are not measured precisely, i.e., 
instead of probability distributions only freqüències of appearances of eigenvalues of these observables are available. 

In the present paper we will focus our attention on the reconstruction of density operators of spin states when mean 
vàlues of a subset of system observables are measured precisely. In this case the estimation of density operators can 
be performed with the help of Jaynes principle of maximum entropy. 



A. Principle of maximum entropy and observation levels 

Let us assume that the state of a spin system, described by the density operator po, is unknown and only expectation 
vàlues G v of observables G u (y = 1, . . . ,n) are available from a measurement. The set of observables is referred to 
as the observation level O []Ï4| . There can be a large number of density operators p which are in agreement with the 
experimental results, i.e., 

Tt{pG v )=G v (u=l,...,n). (2.1) 

If we wish to use only the expectation vàlues G v of the chosen observation level for an estimation (reconstruction) of the 
density op erat or, then we face the problem of selecting one particular density operatorpo out of many p which fulfill 
condition (pTÏ|). To perform this "selection" (i.e., estimation) we note that the density operators under consideration 
do differ by their degree of deviation from pure states. To quantify this deviation an uncertainty measure has to be 



introduced. Following Jaynes |15| one can utilize the von Neumann entropy jïï^ 

S[p] = -Tiip\np). (2.2) 

For pure states 5 = whilc for statistical mixtures of pure states S > 0. 

According to the Jaynes principle of maximum entropy, we have to choose from a set of density operators p which 



fulfill the constraints of Eq.(2.1) the generalized canonical density operator po which maximizes the value of the von 
Neumann entropy^]. In other words, the maximum-entropy principle is the most conservative assignment, in the sense 
that it does not permit one to draw any conclusions unwarranted by the measured data. The generalized canonical 
density operator pç> represents a partially reconstructed (estimated) density operator on the given observation level 
O. The corresponding entropy So = S\po] represents the measure of deviation of the reconstructed state from an 
original pure state. The generalized canonical density operator po takes the form ]Ï4j 



1 Strictly speaking, Wigner functions or wave vectors cannot be measured, they only can be reconstructed from experimental 
data with a help of some inversion procedures. 

2 We note that, in principle, instead of the von Neumann entropy one can utilize another uncertainty measure to distinguish 
between the density operators which fulfill constraints For example, by maximizing the linearized entropy [^| 



n\p] = 1 - Tr(p 



2 v 



one can obtain a partially reconstructed density operator p^ which fulfills conditions (2.1) and simultaneously leads to the 
"maximum mixture", i.e., w[p^] = max. Note that for pure states t] = while for mixtures v > 0. On the complete 
observation level, this "maximum-mixture" principle is equivalent to the maximum-entropy principle but, in general, on the 
reduced (incomplete) observation levels, po / Pq ■ 
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Po = ^- exp (- J2 U , (2.3) 



Zo 
Zo = Tr 



exp 



— 7 \ V G V 



where Zq is the generalized partition function; X v are the Lagrange multipliers which have to be found from the set 



of equations ( |2.l|) 



Any incomplete observation level Oa can be extended to a more complete observation level Ob which includes 
additional observables, i.e., Oa C Ob- Additional expectation vàlues can only increase the amount of available 
information about the state of the system. This procedure is called the extension of the observation level (from Oa 
to Ob) and is usually associated with a decrease of the entropy, as Sb < Sa- We can also consider a reduction 
of the observation level if we decrease number of independent observables which are measured. This reduction is 
accompanied with an increase of the entropy due to the decrease of information available about the state of the 
system. Each incomplete observation level can be considered as a reduction of the complete observation level. In 
what follows we will study a sequence of observation levels in the form 



O a C O b C O c C ... c O comp 
S A > S B > S c > ... > 



(2.4) 



which represent successive extensions of an observation level Oa towards the complete observation level Ocomp- 

Concluding this section we make two remarks. 
(1) Firstly we stress that the reconstruction scheme based on the Jaynes principio of maximum entropy does not 
require any a priori assumption about the purity of reconstructed states, i.e., it can be applied for reconstruction 
of pure states as well as for statistical mixtures. This reconstruction scheme is equivalent to an averaging over the 
generalized g ran d canonical ensemble of all states of the system, under the conditions imposed by the constraints 



given by Eq.(2.1). Within the framework of a geometrical formalism, each state of t he qu antum system is represented 
by a point in the parametric state spacef| Those states which fulfils the constraints (2.1) are represented by a specific 



manifold in the parametric space. From the MaxEnt principle it then follows that the generalized canonical density 
operator is equal to the equally weighted average over all states on this specific manifold. Obviously this average is 
represented by one special point which is associated with the generalized canonical density operator. 
(2) In the case when there is no information available about the preparation of the system, then there is no intrinsic 
way to specify the "minimal" complete observation level. Hcrc by minimal we mcan the complete observation level 
composed of the smallest number of observables. What one can do is to extend systematically observation levels and 
evaluate the von Neumann entropy associated with reconstructed generalized canonical density operators. If at some 
stage of the extension of observation levels the von Neumann entropy becomes zero, it then means that the given 
observation is complete and the pure state of the system is completely "measured". Obviously this does not mcan 
that this observation level is the minimal oneP]. In the case when the measured system is prepared in an unknown 
statistical mixture it is impossible to specify the minimal observation level prior the mcasurement on the complete 
observation level is performed. If this is done then by a sequence of reductions under the condition that the von 
Neumann entropy is unchanged one can specify the minimal observation level. 

In the following sections we will apply the Jaynes principle for the reconstruction of pure spin states. Firstly, 
for illustrative purposes we present the simple example of the reconstruction of states of a single spin- 1/2 system 
with the help of the maximum-entropy principle. Then we will discuss the partial reconstruction of entangled spin 
states. In particular, we will analyze the problem how to identify incomplete observation levels on which the complete 
reconstruction can be performed for the Bell and the Greenberger-Horne-Zeilinger states (i.e., the corresponding 
entropy is equal to zero and the generalized canonical density operator is identical to po). 



III. A SINGLE SPIN-1/2 



Firstly we illustrate the application of the maximum-entropy principle for the partial quantum-state reconstruction 
of single spin-1/2 system. Let us consider an ensemble of spins-1/2 in an unknown pure state \tpo). I n the most 



3 Pure states, which are elements of the generalized microcanonical ensemble, are represented by points on a manifold in this 
state space (such as the Poincare sphere in the case of the spin-1/2). 

4 In fact, in the case of pure states there always exist just one observable (at least in a sense of the Hermitian operator) such 
that the given pure state is an eigenstate of this observable. Unfortunately, it is impossible to specify this operator prior the 
complete reconstruction. 
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general case this unknown state vector \tjjo) can be parameterized as 

IV'o) =cosé»|l) + sin6»|0) (3.1) 

where |0), |1) are eigenstates of the z-component of the spin operator S z = \ò z with eigenvalues —i, i, respectively. 
The corresponding density operator po — \ipo)(ipo\ can be written in the form 

Po = \ (i + n-à) (3.2) 

where I is the unity operator, n = (sin 29 cos <p, sin 29 sin (p, cos 29) ; a — (à x ,à y ,à z ) are the Pauli spin operators which 
in the matrix representation in the basis |0), |1) read 



01 » _ / -i \ „ _ / 1 
1 ) ' ^ ~ l i ' °* ~ l -1 



(3.3) 



To determine completely the unknown state one has to measure three linearly independent (e.g., orthogonal) 
projections of the spin. After the measurement of the expectation value of each observable, a reconstruction of the 
generalized canonical density operator ( |2.3| ) according to the maximum-entropy principle can be performed. In Table 1 

we consider three observation levels defined as = {cr z }, O^g = {à z ,à x } and — {à z ,à x ,à y } = O comp [the 
superscript of the observation levels indicates the number of spins-1/2 under consideration] . 

Using algebraic properties of the <7„-operators, the generalized canonical density operator (2J3) can be expressed as 

Z = 2cosh|A| (3.4) 



1 < 1 
— exp(-A.CT) = - 



cosh |A|7 — sinh|A| j^y 



with A = (X x , X y ,X x ) and |A| 2 = A^ + Xy + X 2 Z . The final form of the po on particular observation levels is given in 
Table 1. The corresponding entropies can be written as 

So = -palnpo - (1 -po)m(l - po) (3.5) 

where po is one eigenvalue of po [the other eigenvalue is equal to (1 — po)] which reads as 



1 + 1(^)1 l + y/(à x ) 2 + (àz) 2 W(^) 2 + K) 2 + (<7,) 2 
Pa = ^ ' Pb = 2 ' Pcomp = 2 ' ^ ' 

It is seen that the entropy So is equal to zero if and only ii po = 1< From here follows that on only the basis 

vect ors |0) and |1) with \{à z )\ — 1 can be fully reconstructed. Nontrivial is Og\ on which a whole set of pure states 
(3.1) with (à y ) — (i.e., <p = 0) can be uniquely determined. For such states Sb — and further measurement of the 
à y on O comp represents redundant (useless) information. 



IV. TWO SPINS-1/2 

Now we assume a system composed of two distinguishable spins-1/2. If we are performing only local measurements 
of observables such as <r^ 1 ' ) ® and ® (here superscripts label the partides) which do not reflect corre- 
lations between the partides then the reconstruction of the density operator reduces to an estimation of individual 
(uncorrelated) spins-1/2, i.e., the reconstruction reduces to the problem discussed in the previous section. For each 
spin-1/2 the reconstruction can be performed separately and the resulting generalized canonical density operator is 
given as a tensor product of particular generalized canonical density operators, i.e., p — p^ ® p^ 2 \ In this case just 
the uncorrelated states |V>o) = | ? /'o 1 ' > ) ® l^o 2 "*) can be fully reconstructed. Nevertheless, the correlated (nonfactorable) 

states \ipo) ^ IV^) ® iW) 2 ^) are °f central interest. 

In general, any density operator of a system composed of two distinguishable spins-1/2 can be represented by a 
4x4 Hermitian matrix and 15 independent numbers are required for its complete determination. It is worth noticing 
that 15 operators (observables) 

{Ò^®ÏV\ï^®G\V\ò^®Gfï} (ii,u = x,y,z) (4.1) 
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together with the identity operator JW ® form an operator àlgebra basis in which any operator can be expressed. 
In this "operator" basis each dcnsity operator can be written as 



j(D 



r(2) 



7(2 



n 



(2) f(i) 



.(2) 



(1) 



r (2) 



(4.2) 



(i) 



ct1 2) ) (/í, v = x,y,z). 



with = (ò% 

Using the maximum-entropy principle we can (partially) reconstruct an unknown density operator po on various 
observation levels. Conceptually the meth od of maximum entropy is rat her straightforward: one has to express the 
generalized canonical density operator (2.3) for two spins-1/2 in the form (|4.2|) from which a set of nonlinear equations 
for Langrange multipliers A„ is obtained. Due to algebraic properties of the operators under the consideration the 
practical realization of this programme can be technically difficult (see Appendix A) . 

In Table 2 we define some nontrivial observation levels. Measured observables which define a particular observation 
level are indicated in Table 2 by bullets (•) while the empty circles (o) indicate unmeasured observables (i.e., these 
observables are not included in the given observation level) for which the maximum-entropy principle "predicts" 
nonzero mean vàlues. This means that the maximum-entropy principle provide us with a nontrivial estimation 
of mean vàlues of unmeasured observables. The generalized canonical density operators which correspond to the 
observation levels considered in Table 2 are presented in Table 3. The signs "©, 0" are used to indicate unmeasured 
observables for which nontrivial information can be obtained with the help of the maximum-entropy principle. 



A. Reconstruction of Bell states 



In what follows we analyze a partial reconstruction of the Bell states (i.e., the most correlated two particle states) 

on observation levels given in Table 2. One of our main tasks will be to find the minimum observation level (i.e., the 

set of system observables) on which the complete reconstruction of these states can be performed. Obviously, if all 

15 observables are measured, then any state of two spins-1/2 can be reconstructed precisely. Nevertheless, due to the 

quantum entanglement between the two partides, measurements of some observables will simply be redundant. To 

find the minimal set of observables which uniquely determine the Bell state one has to perform cithcr a scquence of 

(2) 

reductions of the complete observation level, or a systematic extension of the most trivial observation level O A . 
Let us consider particular examples of Bell states, of the form 

|*(f e ")) = -^[|l,l)+e^|0,0)], p^") = |^(BeH) )( ^(BeH)| ) (43) 

(other Bell states are discussed later). These maximally correlated states have the property that the result of a 
measurement performed on one of the two spins-1/2 uniquely determines the state of the second spin. Therefore, 
these states find their applications in quantum communication systems |Ï8j . In addition, they are suitable for testing 
fundamental principies of quantum mechanics ]Ï9| such as the complementarity principle or local hidden-variable 
theories pp|| . 

(2) 

Let us analyze now a sequence of successive extensions of the observation level O a 

{ A ] C ( g ] C ( ç ] C 0% ] . (4.4) 

(2) 

The observation level O a (see Tab.2) is associated with the measurement of ò z observables of each spin individually, 

(2) 

i.e., it is inscnsitivc with respect to correlations between the spins. On O b both z-spin components of particular 

spins and their correlation have been recorded (simultaneous measurement of these observables is possible because 

(2) 

they commute) . Further extension to the observation level on O c corresponds to a rotation of the Stern-Gerlach 

( 2) 

apparatus such that the x-spin component of the second spin-1/2 is measured. The observation level O d is associated 

with another rotation of the Stern-Gerlach apparatus which would allow us to measure the y-spin component. The 

(2) (2) (2) 

generalized canonical density ope rators on the observation levels O b , O c and O d predict zero mean vàlues for all 



the unmeasured observables (4.1) (see Table 3) 



(2,) 

In general, successive extensions (4.4) of the observation level O a should be accompanicd by a decrease in the 



entropy of the reconstructed state which should reflect increase of our knowledge about the quantum-mechanical 
system under consideration. Nevertheless, we note that there are states for which the entropy remains constant when 
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Og^ is extended towards O^' and O e ' , i.e., the performed measurements are in fact redundant. For instance, this is 



l(2) 



,(2) 



the case for the maximally correlated state (4.3). Here entropies associated with given observation levels read 



S A = 21n2, S B = S c = S D = ln 2, 



(4.5) 



respectively, which mean that these observation levels are not suitable for reconstruction of the Bell states. The 
reason is that the Bell states have no "preferable" direction for each individual spin, i.e., (àp) = for jj, = x,y,z 
and p = 1, 2. 

From the above it follows that, for a nontrivial reconstruction of Bell states, the observables which reflect correlations 
between composite spins also have to be included into the observation level. Therefore let us now discuss the sequence 
of observation levels 



Ol 2) C C O, 



(2) 



l(2) 



(4.6) 



associated with simultaneous measurement of spin components of the two partides [see Table 2] . The corresponding 
generalized canonical density operators are given in Table 3. To answer the question of which states can be completely 

(2) A — r3 

reconstructed on the observation level O e we evaluate the von Neumann entropy ( |2.2| ) of the generalized canonical 
density operator /3g. For the Bell states we find that Se = ~Pe hips — (1 — Pe) hi(l — Pe) where pe = (1 — cos y)/2. 
We can also compare directly pl^ el1 ^ with ç>e- The density operator çr^ in the matrix form can be written as 



píBell) i 



( 1 e~ ü 





\e lv 1 



(4.7) 



(2) 

while the corresponding operator reconstructed on the observation level O e reads 



PE 



( 1 costp 





\ cos ip 1 



(4.8) 



-(2) i "(1) 

al and a x 



f — 0,7T/ 



(2) 



and 



We see that p[^ eU ^ — ps and S[pe] = only if ip — or n which means that the Bell states 
[|1, 1) ± |0, 0)] are completely determined by mean vàlues of two observables a* 

(2) 

that these states can be completely reconstructed on O e . We note that two other maximally correlated states 

1 (2) 

|$±) = ^7= [|0, 1) ± 1 1 , 0)] can also be completely reconstructed on O e . 
i(2) + „ /n( 2 ) 



The extension of O e to O f does not increase the amount of information about the Bell states (fL3) with tp ^ 0,7r. 



(2.) i'i 

For this reason we have to consider further extension of O f to the observation level Oq ' (see Table 2 and Appendix 



-)(2) 



A). In what follows we will show that this is an observation level on which all Bell states (4.3) can be completely 
reconstructed. To see this one has to realize two facts. Firstly, the generalized canonical density operator pq given 
by Eq.(|2.3|) can be expressed as a linear superposition of observables associated with the given observation level, i.e.: 



PG 



\ p,=x,y,z / 
\ fi=x,y,z / 



(4.9) 



where the parameters £ M „ are functions of the Lagrange multipliers X^. Secondly, for Bell states (4.3) the only 



observables which have nonzero expectation vàlues are those associated with Oq\ Namely, {cr)f ® ò 



.(2) 



= 1, 



à x 2 ^) — —(ày 1>( 3óy Z> ) — cos(p and {à x L> ®a 



(!)^^(2) 



ày 1 ^ ®è x 2 ^) — sin<p. It means that all coemcients in the generalized 



canonical density operator pq given by Eq. ( |4.2| ) are uniquely determined by the measurement, i.e., pq = p v . 



(2) 

From the above it follows that Bell states can be completely reconstructed on the observation level O g . On the 

(2) 

other hand, O g is not the mínimum observation level on which these states can be completely reconstructed. The 
minimum set of observables which would allow us to reconstruct Bell states uniquely can be found by a reduction of 
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(2) 

O g . Direct inspection of a finite number of possible reductions reveals that Bell states can be completely reconstructed 
on those observation level which can be obtained from Oq when one of the observables ® (y = x, y, z) is 

(2) 

omitted. As an example, let us consider the observation level O h given in Table 2 which represents a reduction of 
O n when the observable à^®à^ is omitted. Pcrforming the Taylor series expansion of the generalized canonical 
density operator pn defined by Eq. (2.3) one can find that the only new observable ài 1 ^ (g> enters the expression 

for the ph as indicated in Table 3. The coefhcient t in front of <8> can either be found explicitly in a closed 
analytical form (see Appendix A) or can be obtained from the following variational problem. Namely, we remind 
ourselves that the expression (2.3) for pn helps us to identify those unmeasured observables for which the Jaynes 
principle of the maximum entropy "prediets" nonzero mean vàlues. At this stage we still have to find the particular 
value of the parameter t for which the density operator pjj in Table 3 leads to the maximum of the von Neumann 
entropy. To do so we search through the one-dimensional parametric space which is bounded as — 1 < t < 1. To 
be specific, first of all, for t € (—1, 1) we have to exclude those operators which are not true density operators (i.e., 
any such operators which have negative eigenvalues). Then we "pick" up from a physical parametric subspace the 
generalized canonical density operator with the maximum von Neumann entropy. Direct calculation for Bell states 
shows that the physical parametric subspace is reduced to an isolated "point" with t = (ài^ = 1. Therefore we 

conclude that Bell states can completely be reconstructed on Oh ■ Two other minimum observation levels suitable for 
the complete reconstruction of Bell states can be obtained by a reduction of Oq ' when either ® or à\P ® à y 2 ^ 
is omitted. On the other hand, direct inspection shows that a reduction of Oq ^ by exclusion of either aj^ ® à y 2 ^ or 
ày 1 ^ ® àí^ leads to an incomplete observation level with respect to Bell states. 

(2) (2) 

In what follows we discuss briefly two other observation levels 0\ and Oj which are defined in Table 2. The 

observation level O \ serves as an example when one can find an analytical expression for the Taylor series expansion 
of the canonical density operator pi ( |2.3| ) in the form_([T^). The coefficients (functions of the original Lagrange 
multipliers) in front of particular observables in Eq.(4.2) can be identified and are given in Table 3. Problems do 

(2) (2) i i 

appear when 0\ is extended towards O T . In this case we cannot simplify the exponential expression (2.2) for pj 



and rewrite it analytically in the form ( |4.2| ) as a linear combination of the observables ( J4. l| ) . In this situation one 
should apply the following procedure: firstly, by pcrforming the Taylor-series expansion of the pj to the lowest orders 
one can identify the observables with nonzero coefficients in the form (4.2). Namely, for pj the additional observables 
àz <8> &x , óí 1 ^ ® 5* and Oy ® ày 2 ^ appear in addition to those which form [see Table 3]. The corresponding 
coefficients u,v,w G (—1,1) form a bounded three-dimensional parametric space (u,v,w). In the second step one 
can use constructively the maximum-entropy principle to choose within this parametric space the density operator 
with the maximum von Neumann entropy. The bàsic procedure is to scan the whole three-dimensional parametric 
space. At the beginning, one has to select out those density operators (i.e., those parameters u,v,w) which posses 
negative eigenvalues and do not represent genuine density operators. Finally, from a remaining set of "physical" 
density operators which are semi-positively defined the canonical density operator pj with maximum von Neumann 

(2) (2) 

entropy has to be chosen. For a completeness, let us notice that for Bell states the observation levels O) and O ) 

(2) 

are equivalent to O e , i.e., pi — p,j — Pe- 

(2) 

In this section we have found the minimum observation levels [e.g., O h } which are suitable for the complete 
reconstruction of Bell states. These observation levels are associated with the measurement of two-spin correlations 
&X (g) ctÍ 2 -*, òy^ ® <ri 2 · ) and two of the observables &u <8 (y = x, y, z). Once this problem has been solved, it is 
interesting then to find a minimum set of observables suitable for a complete reconstruction of maximally correlated 
spin states systems consisting of more than two spins-1/2. In the following section we will investigate the (partial) 
reconstruction of Greenberger-Horne-Zeilinger states of three spins-1/2 on various observation levels. 



V. THREE SPINS-1/2 



Even though the Jaynes principle of maximum entropy provides us with general instructions on how to reconstruct 
density operators of quantum-mechanical systems practical applications of this reconstruction scheme may face serious 
difficulties. In many cases the reconstruction scheme fails due to insurmountable technical problems (e.g. the system 
of equations for Lagrange multipliers cannot be solved explicitly) . We have illustrated these problems in the previous 
section when we have discussed the reconstruction of a density operator of two spins-1/2. Obviously, the general 
problem of reconstruction of density operators describing a system composed of three spins-1/2 is much more difficult. 
Nevertheless a (partial) reconstruction of some states of this system can be performed. In particular, in this section 
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we will discuss a reconstruction of the maximally correlated three spin-1/2 states - the so-called Greenberger-Horne- 
Zeilinger (GHZ) state @ : 



|*^) = -^[|l,l,l) + e i ^|0,0,0)J , pf HZ ^ = \^ ¥ )^ 



(5.1) 



Our main task will be to identify, with the hclp of the Jaynes principle of maximum entropy, the minimum observation 
level on which the GHZ state can be completely reconstructed. 

We start with a relatively simple observation level O b such that only íwo-particle correlations of the neighboring 
spins are measured, i.e. 



Of = {a« ® af ® /( 3 ),/« © fi-f © ai 3 )}. 
The generalized density operator associated with this observation level reads 

/5 B = i [l« ® j< 2 > ® /i 3 ) + (ài 1 ) © àf © Í&)àW g, 42) g, ; ( 3) 

© (a« ® ® /( 3 ))(/W <rf cfV^ 7 2 ai 3 )" . 



(5.2) 



(5.3) 



where '©" indicates a prediction for the unmeasured observable. In particular, for the GHZ states (5.1) we obtain 
the following generalized canonical density operator 



JGHZ) 

Pb 



1 r 



f(l) ffl f( 2 ) 



f(3) 



+ ai 1 ) ai 2 ) © + I« ai 2 ) a^ © <7« © J 2 
= i|l,l,l)(l,l,l| + i|0,0,0)<0,0,0|. 



(5.4) 



The reconstructed density operator ff^ HZ ^ describes a mixture of three-particle states and it does not contain any 

information about the three-particle correlations associated with the GHZ states. In other words, on Og the phase 
information which plays essential role for a description of quantum entanglement cannot be reconstructed. This is 
due to the fact that the density operator p^ HZ ^ is equal to the phase-averaged GHZ density operator, i.e. 



~(GHZ) 

Pb 



1 

2^ 



pf HZ) &p. 



(5.5) 



Because of this loss of information, the vonNeumann entropy of the state p^ HZ ^ is equal to ln 2. We note, that when 
the GHZ states are reconstructed on the observation levels O^} — {à^ © àffi © j( 3 ), JÍ 1 ) © © a^} (p = x,y), 



then the corresponding reconstructed operators are again given by Eq.( p.4| ). These examplcs illustrate the fact that 
three-particle correlation cannot be in general reconstructed via the measurement of two-particle correlations. 

To find the observation level on which the complete reconstruction of the GHZ states can be performed we recali 
the observables which may have nonzero mean vàlues for these states. Using abbreviations 



= (4 X) ® ® Í W ), $H»3 = ® ® àí 3) ), ^ = (fi-W © ï {2) © àl^), 



we find the nonzero mean vàlues to be 



U 1 ) 



-(2) 



»(3)v 



(|U,i/,w = x,y,z), 



(5.6) 



£ziZ2 = ^2 2 2 3 = ÇZ1Z3 = 1) 

CyiV2X3 CxiJ/2J/3 CyiX2V3 
ÇxiX2X3 COS y?. 

Cyiy2V3 = —suiip. 



sm (/?, 

— COSí/3, 



(5.7) 
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We see that for arbitrary ip there exist non-vanishing three-particle correlations Cah^^s- The observation level which 
consists of all the observables with nonzero mean vàlues is the complete observation level with respect to the GHZ 
states. Our task now is to reduce this set of observables to a minimum observation level on which the GHZ states 
can still be uniquely determined. In practice it means that each observation level which is suitable for the detection 
of the existing coherence and correlations should incorporate some of the observables with nonzero mean vàlues. The 
other observables of these observation levels should result as a consequence of m utual tensor products which appear 
in the Taylor series expansion of the generalized canonical density operator ( |2.3| ). It can be seen by direct inspection 
of the finite number of possible reductions that the minimum set of the observables which matches these requirements 
consists of two two-spin observables and two three-spin observables. For the illustration we consider the observation 
level 



J(3) J(D 



r (2) 



f(3) 



><1) 

y 



r(3) 



}• 



(5.8) 



In this case the exponent C of the generalized canonical density operator pc = exp(— C)/Zc [see Eq.(2.3)] can be 



rewritten as C = C\ + C2 with C\ — "fuàl 



pày L> ® ày ® à y 0> . The operators C\, C2 commute and further calculations are straightforward. After some àlgebra 
the generalized density operator pc can be found in the form 

1 



;.(3) 



ctÍ 2) ® /( 3 ) + <8 ò z 2) ® àf ] and C 2 



,(1) 



Pc = - ® /( 2 ) ® m + ç zlZ2 &p ® &w ® m + 6 2 , 3 Í (1) ® à? ] ® <^ 3) 

o L 



(5.9) 



r (3) 



4(1) 



r(2) « A(3) 



For the GHZ states the vonNeumann entropy of the generalized canonical density operator pc is equal to zero, from 



which it follows that pc — p^ HZ ^ [see Eq.(S.l)], i.e., the GHZ states can be completely reconstructed on O, 



c 



Moreover, the observation level Oc represents the minimum set of observables for complete determination of the 
GHZ states. 



VI. CONCLUSIONS 



We have investigated the problem of a (partial) reconstruction of correlated spin states on different observation 
levels. We have found the minimal set of observables for the complete reconstruction of the most correlated states for 
systems composed of two and three spins-1/2, i.e., Bell states and GHZ states. Direct generalization to systems of 
more spins-1/2 is possible. 

The concept of observation levels and the maximum-entropy principle is a powerful tool which can be used also for 
other physical system, e.g., for the reconstruction of the states of a monochromatic light-field ]2Ï] ]. We recali that this 
reconstruction scheme is based on the knowledge of the exact mean vàlues of given observables or their probability 
distributions (see Appendix A) . Theoretically, this means that an infinite number of measurement over an ensemble of 
identically prepared system has to be performed in order to obtain those mean vàlues which are needed. In practice, 
if the number of measurements is sufficiently high, then the mean vàlues can be considered to be measured precisely 
enough and the Jaynes principle can be applicd for a state reconstruction. On the other hand, if just few measurements 
are performed, then the mean vàlues of the considered observables are not known and the Jaynes principle cannot be 
used. 

In this case another reconstruction scheme has to be applied. In particular, in the case of a small number of 
measurements the Bayesian reconstruction scheme can be effectively utilized. We will address the problem of a 
reconstruction of correlated spin systems based on Bayesian methods elsewhere p3| . 
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APPENDIX A 



Conceptually the reconstruction scheme based on the Jaynes principle of the maximum entropy is very simple. 
On the other hand particular analytical calculations can be difíicult and in many cases cannot be performed. In 
this appendix we present explicit calculations of generalized canonical density operators (GCDO) and corresponding 

(2) (2) 

entropies for two observation levels O g and O h defined in Table 2. 



(2) 

A. 1. Observation level O a 

Let as assume the observation level O g given by the set of observables {à z 
à^; ày 1 ^ ® à y 2 ^}. In this case the GCDO reads 

p G = ±- exp (-É) (A.l) 

where 



Z G = Tr 



exp 



(-È)] (A.2) 



is the partition function. Here we have used the abbreviation 

È = x zz ài 1} ® ^ 2) + K x ài 1] ® ^ + \ xy ò^ ® + \ yx d y v ® + Xyyà^ ® . (A.3) 
The corresponding entropy has the form 

Sg = hlZ<3 + ^zz£,zz + ^xx£,xx + ^xy^xy + ^yxÇyx + -^yy£,yy> (A-4) 

Using the algebraic properties of the operators associated with the given observation level we find the GCDO (A.l) 
to read 



® -< 2) + e B x^ J ® 4 2) + ® ^ 2) 



(A.5) 



where we use the notation 



tu» = ® àl 2) ) , (n,v = x,y,z). (A.6) 



Now we express the entropy as a function of expectation vàlues of operators associated with the observation level 

(2) (2) (2) (2) 

O g . With the hclp of this entropy function we can perform reductions of O g to the observation levels O h , O f 
and Og\ In order to perform this reduction we express A 

(xv in Ecj. (A.4) as functions of the expectation vàlues 

To do so we utilizc the rclation 

- - d -B- ™ 

The partition function Zq can be found when we rewrite the operator E in Eq.(A.4) as a 4x4 matrix: 

(a d*\ 
s 7 -i o • (-> 
d a J 

where we used the abbreviations 
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The powers of the operator É can be written as 



E n = 



with the matrix elements given by the relations 













\ 





^22 


-^23 










-^32 


-^33 







Ui? 








^44 


/ 



= BW = I[(a+|d|) n + (a-|d|n 



(n) 
14 
(n) 
22 
(n) 
23 
(n) 
32 
.(n) 
41 



i[(«+Mir-( a -Min F . 



£&M[(-«+i&i) n +(-«-i&ir 

i[(- + |6|) n -(-o- 



(n)* 



Now we find 



exp 



(-*) 



e a cosh 



V- e - a sinh(|d|)^ 



e a cosh |6| 
-e«sinh(|6|)^ 



\d\ 





-e a sinh(|6|)^ 
e a cosh |6| 




- e - a sinh(|d|)^^ 



e- a cosh|d| j 



from which we obtain the expression for the partition function Zq 

Z G = 2e- a cosh |d| + 2e a cosh |ò| . 

For the expectation vàlues given by Eq.(A.7) we obtain 

Ç zz = JL [2e' a cosh \d\ - 2e a cosh |6|] ; 

2e- a sinhfldQ |ij (X xx - X yy ) + 2e a sinh(|6|)^ (X xx + X yy ) 
2e- a smh(\d\) |ij (X xy + X yx ) + 2e a sinh(|ò|)^ {X xy - X yx ) 
2e- a sinhfldQ^ {X xy + X yx ) - 2e a sinhflòQ^ (X xy - X yx ) 
-2e- a sinh(|d|)^ T {X xx - X yy ) + 2e a sinh(|ò|)^ (X xx + X yy ) 



i 

Z G 
1 

Z G 



p - I 



Z G 

P = L· 

Çyy z G 

If we introduce the abbreviations 

B ^ xx -+- £ >yy i {^ xy 
thcn with the hclp of Eq. (A.14) we obtain 



S,yx) ) D — ^ XX l^yy + Í (Ç Xy + (.yx) 



Taking into account that 



we hnd 



\B\ = e a sinh(|6|), \D\ 

Zg 



B b D 

lbï""W' W\ 



-^e-«sinh(|,|)^. 



— e - a sinh(M|) 



Now we introduce four new parameters Mj 

Aíi = l+^„ + |D| , M 2 



1 + 6 



\D\ 
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M 3 = 1-Ç ZZ + \B\ , M A = 1-Ç ZZ -\B\ (A.19) 

in tcrms of which we can express the vonNeumann entropy on the given observation level. Using Eqs. (A.13), (A.14) 
and (A.17) we obtain 



M i = ^ ex P(-«+MI) - M 2 = ^exp(-a-|d|) , 
M 3 = ^exp( a +|ò|) , M 4 = £exp(a-|6|). lA ' 2Uj 

The Lagrange multipliers A/^ can be expressed as functions of the expectation vàlues 



, s /M 3 M 4 \ 4 „ 1N fM 3 \ 2 „ lx /MA 2 

Aftcr inscrting these expressions into Eq. (A.13) we obtain for the partition function 

Z G = r- (A.22) 

(M1M2M3M4) 3 



When we insert Eqs. (A.18), (A.21) and (A.22) into Eqs. (A.l), (A.4) and (A.12) then we find both the entropy 
and the GCDO 



4 V 4 



PG 







+ Czz 








D* 


1 





1-6, 


B* 





4 







B 


l-ízz 







V 


D 






i + e 



(A.24) 



as functions of the expectation vàlues Finally, wc can rewrite the reconstructed density operator (A.24) in terms 
of the spin operators (see Tab.3). 



dS G 1 (M X M 2 \ 
Kz = 1 ^ ± = - 7 In 7/7/ = (A.25) 



A. 2. Observation level 

(2) (2) 

The GCDO on the 0%' can be obtained as a result of a reduction of the observation level O g ' . The difference 

between these two observation levels is that the does not contain the operator ®ò z 2 \ i.e., the corresponding 
mcan valuc is unknown from the measurement. 

(2) (2) 
According to the maximum-entropy principle, the observation level O h can be obtained from O g by setting the 

Lagrange multiplier X zz equal to zero. With the help of the relation [see Eq.(A.7)] 

dS G _ 1 / M\M; 

we obtain 

MiM 2 = M3M4. (A.26) 
From this cquation we find the "predicted" mean valuc of the operator ài^ <g> à z 2 ^ (i.e., the parameter t in Table 3) 

^ = \{\D\ 2 -\B\ 2 )=t. (A.27) 
Taking into account that the parameters \B\ and \D\ read 

\B\ = (Çxx + £yy) + (£xy — £yx) j \D\ = (Çxx ~ Çyy) + (Çxy + ' (A.28) 
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we can express the predicted mean value Ç zz as a function of the measured mean vàlues £ xx , £ xy , Ç yx and Ç yy : 

Çzz — (S,xy Çyx - Íxx^yy) • (A.29) 

When we insert Eq. (A.27) into Eq. (A.19) we obtain: 

Mi = N X N 2 , M 2 = N 3 N 4 , M . = A1A3 , M 4 = N 2 N 4 , (A.30) 
where the parameters Ni are defined as 

N 1 = 1 + ^(\D\ + \B\) , N 2 = l + ±(\D\-\B\) , 

N 3 = N 4 = í-l(\D\ + \B\) . (A.31) 

In addition, from Eqs. (A.30) and (A.23) we obtain the expression for the vonNeumann entropy of the density 

(2) 

operator reconstructed on the observation level O h . 

Finally, from Eqs. (A.28) and (A.24) we find the expression for the GCDO on the observation level Ojp (see Tablc 
3): 



PH = \[ï {1) ®ï {2) + (t, xy íyx~íxxíyy)àz 1) ®àz 2) 



+ ÇxxàP <g> Ò x 2) + Ç x yà x ^ <g> (Ty 2) + Çyxày^ ® Òx^ + ÇyyÒy^ ® Uy^] . 
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reconstructed density operator 
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Pa = \ [l + n z & z ^j 


^B 
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PB = \ [i + n z ír z + n x a x ^j 


o (1) 

^comp 
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• 


• 


Pcomp = \ [ï + n z ír z + n x ír x + ny&y^j 



Table 1 

sociated wi 

indicate which observables constitute a given observation level. We also present explícit expressions for the reconstructed density operators px', Pb and p CO mp 



In this table we present three observation levels ; \ and Ocoïnp associated with a measurement of the particular spin-1/2 operators. Bullets (•) in the table 
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Table 2. 

We present a set of observation levels on which the density operators of two spins-1/2 can be partially reconstructed. Bullets (•) in the table indicate which observables 
constitute a given observation level while empty circles (o) denote unmeasured observables (i.e., these observables are not included in the given observation level) for which 
the maximum-entropy principle "predicts" nonzero mean vàlues. 
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reconstructed densíty operator 


^ K 


p A = \ (/(i)/(2) + nP&PlW + nPlW&P 8 nPnP&P&P) = \ + n { P (/< 2 ) + 4 2) 4 2) ) 


Cl • ' 


p B = \ (/ (1) / (2) + 4 ví^/í 2 ) + nr , /( i )4 J +6,4 vr^j 


Cl • ' 


Pc = \ [i^lW + 4 V^/í 2 ) + nY'lWvr + ni A) I^4 A> + 6,4 + 6,4 V^j 


Cl • ' 


PD = 4 (/(^/W +4 'V^/í 2 ) + 4 J / (1) 4 ' + 4 ^'I^trr +4 ^ ( '4 +6,4 ><?¥> + 6,4 J 4^ +^ y 4 1 V^J 




Pe = | (/ (1) / (2) + 6,4 1) 4 2) +U*Í ) *Í ) 6 6,6,4 '4 




/5 F = \ [i^iw +6,4 } 4 } +6,4 } 4 } +6^4 4 ') 


^(2) 


pg = \ (/ (1) / (2) + 6,4 J 4 J +6,4 J 4 J + 6^4 4 +6 y 4 J 4 +6,4 4 J j 


o (2) 


í _ 1 ( r(l) r(2) , í ^.(1)^.(2) .í ^.(1)^(2) ,t ^M 2 ) , ^UU(2) ffi ,-(l)-(2)\ . f-f f f f 
kH — 4 l ^ íri^x u r T syy u y u y ^ Çxy u x u y \ Çyx u y u x W z u z J j ^ — ÇxyÇyx ÇxxÇyy 


(2) 






P3 = \ (/(D/( 2 ) + 4^4^/í 2 ) + nPírWiW + n^/í 1 )^ + n^I^ + 6,4 1) 4 2 > + 6,4 1) 4 2) 8 8 „*W*< 2 > 6 v,»™»™) 



Table 3 

We present explícit expressions for the reconstructed density operators px of two spins-1/2 on the observation levels denoted in Table 2. We use the notation 
= (4 P ' 1 ) (a* = Z ) X )V') P = 1)2) and 6^ = (4 1 ' ® 4 2 ') with fí,v = z,x,y. The signs © and are used to indicate unmeasured observables for which nontrivial 
information can be obtained with the help of the maximum-entropy principle. 
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